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O , Abstract 

^ ' In this correspondence our aim is to use some tight lower and upper bounds for the differential 
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H 



quaternary phase shift keying transmission bit error rate in order to deduce accurate approximations 
for the bit error rate by improving the known results in the literature. The computation of our new 
approximate expressions are significantly simpler than that of the exact expression. 
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(N ■ I. Introduction 

en . _ 

O ■ Recently, Ferrari and Corazza [T] have studied the performance analysis of the differential 

(N 

quaternary phase shift keying (DQPSK) transmission with Gray coding over the additive white 

Gaussian noise (AWGN) channel. We note that the bit error rate (BER) in the case of Gray 

^ ■ 
;h i coding can be written in terms of the Marcum Q-function and the modified Bessel function of 

c^ ■ 

the first kind and zero order. In [1] the authors used some bounds for the Marcum Q-function in 

order to propose some simple, but accurate approximations for the BER. By using some other 

bounds for the Marcum Q-function (deduced by Wang and Wu [|3|, and Baricz and Sun [2J), 
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very recently Sun et al. 0] proposed another approximation for BER. In this correspondence we 
make a contribution to the subject by using some new results on Marcum Q-function: we deduce 
a new tight upper bound for the BER and combining this with some existing tight bounds we 
construct some new accurate approximations. These results improve and complement the results 
of Ferrari and Corazza Hlj, and also of Sun et al. [j4]. In Section |Il] we discuss the bounds for 
BER of DQPSK, in Section [nil we present five new approximations for BER and we compare 
them with the similar ones deduced from inequalities for the Marcum Q-function. Finally, in 
Section |W] the conclusion of the paper is given. 

II. Bounds for BER of DQPSK 

It is known that for the DQPSK transmission with Gray coding over an AWGN channel the 
BER can expressed as follows [21: 

BER ^g(a, 6) -^/o(a6)e-^, (1) 



where a = y 7(2 — a/2), b = y 7(2 + y/2) and 7 is the bit signal-to-noise ratio (SNR). Here 
Q{a,b) stands for the Marcum Q-function, defined by 



fOO 

A 



Q{a,b) = xe 2 lQ(ax)dx, 
Jb 

where 6 > 0, a > and Iq is the modified Bessel function of the first kind and zero order. Note 
that the Marcum Q-function play an important role in the studies of digital communications over 
fading channels [^. 

Now, let us consider the complementary error function, defined by 



2 



erfcfx) = — ^ / e ^^dt, 

vr 



00 
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and for shortness let us introduce the notation 



00 



e(a, 0) = erfc — ^ = —= / e at. 
V V2 / V^^^ 

Recently, Ferrari and Corazza [[D by using their bounds from [[71 (see also [_8J) for the Marcum 

Q-function _ 

TT 6/0 (afe) 

2~e" 



Qia,b)>J-^^eia,b), (2) 






+ a\l-e{a,b) 



(3) 
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where 6 > a > 0, deduced the following bounds for BER of DQPSK with Gray coding: 



BER > Io{ab) 



BER < Io{ab) 



TT a . 1 _<£+&£ 



u, 



Ui. 



It is important to mention here that motivated by the above results of Ferrari and Corazza, very 
recently Wang and Wu [3] proved that 

/vT hIo{ah) 



Q{a,b) > 

while Baricz and Sun flU proved that 

Io{ab) 



2e 



ab 



-ab 



E{a,b), 



Q{a,b)< 



_{b-af 

e 2 



pab _|_ p—ab 

where in both of the inequalities 6 > a > and 






b+a 
V2 



E{a,b) = e{a,b) — e{—a,b) = —= / e dt. 

Now, by using the above lower and upper bounds for the Marcum Q-function, Sun et al. |J4| 
deduced the following: 



BER > Io{ab) 



BER < Io{ab) 



b 1 _s^+iP 
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We note that the lower and upper bounds for the Marcum Q-function of Wang and Wu S and 
of Baricz and Sun [|2l, mentioned above, are tighter than ^ and ([3]), as it was pointed out in 
[El, S. This in turn implies that the bounds L2 and U2 are tighter than the bounds Li and Ui, 
that is, we have L2 > Li and U2 < Ui. 

Now, let us consider the following inequality 



Q{a,b] < —r — 



e(a,6) + (e'^'' + Ao)e — ^ 



which holds for all 6 > a > and it was deduced by Baricz [|9l. Here Aq = e''''(/o(po)/A(Po) ^ 
1) ~ 3.03442206626763, po ^ 1.54512596391949 is the unique simple positive root of the 
equation {x + l)/i(a;) = xIo{x), and h stands for the modified Bessel function of order 1. We 
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note that in the above inequality the constant Aq is best possible, that is, cannot be replaced by 
any larger constant. Since the above inequality is an improvement of Q, its equivalent form 



BER < Io{ab) 



-K a , ,, 1 
r ;— e(a, 0) H — e 



U. 



clearly provides a better upper bound than (|3}, that is, we have U3 < Ui. The lower and upper 
bounds Li, L2, Ui, U2 and U3 are shown in Fig. [T] together with the exact BER as functions of 
the SNR on the interval (0, 1.5). Surprisingly, all of these bounds are extremely tight. 



III. Approximations for BER of DQPSK 

Observe that there is an interesting formal symmetry between the bounds Li and Ui, as well 
as between the bounds L2 and U2. Based on the tight bounds Li and Ui, and in view of the 
formal symmetry, Ferrari and Corazza [1] pointed out that it is natural to approximate the BER 
simply by considering the arithmetic mean of the quantities Li and Ui . In this spirit, they derived 
the following approximate expression for BER: 

BERi =^(Li,Ui;0.5), 



I.e. 



BERi 



3 ^^o(a&)e(a,6), 



where A{x,y]uj) = wx + (1 — ijj)y is the weighted arithmetic mean of x and y with weights 
cj, 1 — a; > 0. By using the same idea, it is natural to propose the following new approximate 
expressions for BER: 

BER2=^(L2,U2;0.5) 



and 



that is. 



and 



BER3^^(L2,U3;0.5) 
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We note that these approximations are better than BERi. A similar approximation to that of 
BER2 was proposed recently by Sun et al. [lU in the form 



e 2 1 

BER4 = , + 



, , , , w ? + \/- 1 E{a,b), 

VS^ 4 \ V & V a ' ^ ^ 



however, this is better than BERi only for 7 G (0.9,1.4). Motivated by the above results, in 
this correspondence our aim is to propose the following approximations based on better weight 
functions cjj, i G {5, 6, 7} : 

BER5^^(Li,Ui;a;5), 

BER6^^(L2,U2;a;6), 

BER7^^(L2,U3;a;7), 



where 



^5(7) 



0.65^ 



0.5 + 1.1- 



= 277 



if 7 < 1 
h ifl<7 



^6(7) 



1 



e"^ -0.25 + 0.5, ifO<7<l 

f -1.15 + 0.5, ifl<7<5 



(7+0.5)1-5 

^(1 + 7)^^-0.65 + 0.5, if5<7 



(1-7)2-0.95, ifO<7<l 

^7(7) = I 0.5 - 1.4 - e"^'' + 0.02, if 1 < 7 < 8 • 

^ + 0.5, if 8 < 7 
These approximations can be rewritten as follows 

BER5 = a;5(7)Li + (1 - W5(7))Ui, 

BERe = a;6(7)L2 + (1 - W6(7))U2, 

BER7 = a;7(7)L2 + (1 - W7(7))U3. 

It is important to note here that the computations of approximate BER expressions BERj, where 
i G {1,2,. ..,7}, completely avoid the computation of the Marcum Q-function. Tables H and HI] 
contain an estimate of BER by using the Matlab function marcumq (BER) and the estimates 
based upon BER,;, where i G {1, . . . , 7}. As we can see the approximations BERj, where i G 
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{5,6,7}, are better than the approximations BERj, where i G {1,2,3,4}. Moreover, the novel 
approximate expressions BERj, where i E {5, 6, 7}, are consistent with the results of Weinberg 
IfTOl . obtained by using of Monte-Carlo estimators. We note that in our further calculations we 
considered the values obtained by the marcumq function as the "exact" values for BER. 

We give a hint about the construction of the weight functions. By using the values of BER, 
Li, L2, Ui, U2 and U3 first we calculated numerically the values of the weight functions. The 
representation of the initial numerical values for cuq can be seen in Fig. [21 This looks like the 
probability density function of the normal distribution but having a heavy tail, so we started 
to search for expression involved in probability density functions with heavy tails (eg. Student, 
Levy, a— stable). By using the same idea the parameters in the expression of the functions coi, 
i E {5, 6, 7}, were obtained by numerical experimentations. 

Now, in order to estimate the tightness of the new approximate expression and to compare 
with the tightness of the known approximate expressions, for i E {5, 6, 7} let us consider the 
relative errors Si = (BERj — BER) /BER. Table [ni] contains the relative errors £j, where i E 
{5, 6, 7}. Finally, in Fig. [3] and |4] the exact BER and the approximate expressions BERj, where 
i E {1, . . . , 7}, are shown as functions of the SNR on the intervals (0, 2] and [2, 4]. We find that 
our new approximate expressions are quite accurate for the whole region of SNR. 

IV. Conclusion 

In this correspondence we deduced a new tight upper bound for the BER of DQPSK with 
Gray coding over the AWGN channel by using a new tight bound deduced for the Marcum 
Q-function. This and the bounds from [J2l show that approximations of BER for 7 = 1 obtained 
by using adaptive Simpson quadrature (see [lOJ) and some previous approximations are not in 
the interval determined by the lower and upper bounds. Although the estimations from [H] and 
JH are the best possible estimations for large values of SNR, by using constants in the convex 
combinations of the lower and upper bounds, they can be improved by considering weight 
functions in the calculation of approximations. The advantage of this technique is that we can 
obtain better approximations also for small values of SNR. Based on numerical experiments we 
proposed five new approximate BER expressions for the AWGN channel, which use the best 
known bounds and which are accurate in the whole region of SNR. 
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Fig. 1. BER with lower and upper bounds 




Fig. 2. The graph of the weight function ojq 
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Fig. 3. Approximations for BER on [0, 2] 
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Fig. 4. Approximations for BER on [2, 4] 
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10 



7(dB) 


BER 


BERi 


BER2 


BERs 


1 


1.639 X 10"^ 


1.739 X 10"* 


1.731 X 10"* 


1.556 X 10"* 


2 


7.161 X 10"^ 


7.324 X 10"* 


7.322 X 10"2 


7.007 X 10"^ 


3 


3.422 X 10"^ 


3.458 X 10-2 


3.458 X 10"2 


3.416 X 10"^ 


4 


1.701 X 10"^ 


1.711 X 10"2 


1.711 X 10"2 


1.706 X 10"2 


5 


8.648 X 10"'^ 


8.683 X 10"* 


8.683 X 10"* 


8.677 X 10"* 


6 


4.461 X 10"^ 


4.474 X 10"* 


4.4745 X 10"* 


4.473 X 10"* 


7 


2.325 X 10"^ 


2.330 X 10"* 


2.3308 X 10"* 


2.33072 X 10"* 


8 


1.221 X 10"^ 


1.224 X 10"* 


1.22405 X 10"* 


1.22404 X 10"* 


9 


6.459 X 10-* 


6.468 X 10"* 


6.46883 X 10"* 


6.46881 X lO"* 


10 


3.431 X 10""* 


3.435 X lO""* 


3.43588 X 10""* 


3.43588 X 10""* 


11 


1.830 X 10""* 


1.832 X lO""* 


1.83249 X 10""* 


1.83249 X 10"* 


12 


9.798 X lO"-""' 


9.807 X 10"-' 


9.80723 X lO"-""' 


9.80723 X 10"^ 



TABLE I 
BER AND THE ESTIMATES BERi, BER2 AND BER3 
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11 



7(dB) 


BER4 


BER5 


BERe 


BER7 


1 


1.484 X 10"^ 


1.677 X 10"^ 


1.6383 X 10"^ 


1.645 X 10"^ 


2 


6.908 X 10"^ 


7.179 X 10"2 


7.162 X 10"2 


7.133 X 10"2 


3 


3.348 X 10-2 


3.423 X 10-2 


3.4226 X 10-2 


3.4233 X 10"2 


4 


1.677 X 10^2 


1.7014 X 10"2 


1.7017 X 10"2 


1.7036 X 10-2 


5 


8.5931 X 10"^ 


8.6500 X 10"^ 


8.6500 X 10"'^ 


8.6593 X 10"^ 


6 


4.4788 X 10"^ 


4.4624 X 10"^ 


4.4616 X 10"^ 


4.4651 X 10"^ 


7 


2.3741 X 10"^ 


2.3262 X 10"^ 


2.3257 X 10"^ 


2.3267 X 10"^ 


8 


1.2841 X 10"^ 


1.2222 X 10"^ 


1.2219 X 10"'^ 


1.2218 X 10"'^ 


9 


7.1447 X 10"" 


6.4613 X 10"" 


6.4597 X 10-" 


6.4594 X 10"" 


10 


4.1463 X 10"** 


3.4326 X 10"" 


3.4318 X 10"" 


3.4317 X 10"" 


11 


2.5591 X 10"" 


1.8311 X 10"" 


1.8307 X 10"" 


1.8306 X 10"" 


12 


1.7151 X 10"" 


9.8011 X 10"^ 


9.7990 X lO"-""' 


9.7989 X 10"^ 



TABLE II 
Estimates of BER by using BER4, BER5, BERe and BER7 
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12 



7(dB) 


£5 


£6 


£7 


1 


2.31 X 10^2 


-4.51 X 10"* 


3.86 X 10"^ 


2 


2.52 X 10"^ 


1.96 X 10"* 


-3.79 X 10"-' 


3 


1.17 X 10""* 


-1.36 X 10"'^ 


1.94 X 10"* 


4 


9.01 X 10"^ 


2.66 X 10"* 


1.38 X 10"^ 


5 


1.96 X 10"* 


1.93 X 10"* 


1.26 X 10"^ 


6 


2.52 X 10""* 


8.90 X lO""' 


8.65 X 10"* 


7 


2.72 X 10"* 


4.55 X lO"'' 


5.04 X 10"* 


8 


2.72 X 10""* 


2.54 X 10"''' 


-6.42 X 10"^ 


9 


2.62 X 10"* 


1.53 X lO"'' 


-3.32 X 10"^ 


10 


2.48 X 10"* 


1.00 X 10"''' 


-1.66 X lO"-""' 


11 


2.32 X 10"* 


7.13 X 10"" 


-6.71 X 10"" 


12 


2.17 X 10"* 


5.44 X 10"" 


-3.87 X 10"^ 



TABLE III 

Relative ERRORS £i ofthe approximate expressions BERi, where i e {5,6,7}. 
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